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Abstract 

The response of an infinite, periodic, insulating, solid to an infinitesimally 
small electric field is investigated in the framework of Density Functional The- 
ory. We find that the applied perturbing potential is not a unique functional of 
the periodic density change : it depends also on the change in the macroscopic 
polarization. Moreover, the dependence of the exchange-correlation energy on 
polarization induces an exchange-correlation electric field. These effects are 
exhibited for a model semiconductor. We also show that the scissor-operator 

technique is an approximate way of bypassing this polarization dependence. 
PACS numbers: 71.10.-hx; 77.20. -Fy; 78.20.Bh 
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Density-Functional Theory (DFT) has become the standard method for first- 

principles calculations of the properties of solids. Within the Local Density Approxima- 
tion (LDA) 0, the accuracy obtained for a wide set of properties is of the order of a few 
percent 0. Notable exceptions to this gratifying picture are the cohesive energy of solids, 
the energy gap of semiconductors, and their dielectric susceptibility. The first of these fail- 
ures is attributed to the LDA [^, while the second comes from a discontinuity in the DFT 
exchange-correlation potential when an electron is added across the gap @]. 

The third failure is striking in that it is not easy to see whether it comes from the 
LDA approximation, or from some fundamental feature of DFT. At first sight, the dielectric 
susceptibility, being the second derivative of the ground state energy with respect to an 
applied electric field, should be obtained exactly within DFT. Although the corresponding 
potential change is linear in space, and breaks the periodicity of the solid, techniques were 
found to bypass this problem Unfortunately, the LDA results [^j-^ are usually much 

larger (more than 10%) than the experimental data. Dal Corso, Baroni, and Resta ^ 
argued that any attempt to improve the LDA calculated value should go in the direction 
of better approximations to the true functional. The latter authors tried to compute the 
dielectric response of Si in the Generalized Gradient Approximation (GGA) but with 
limited success. 

Earlier, Levine and Allan had proposed a simple "LDA + scissor correction" approach to 
the dielectric tensor , in which a constant shift A is imposed on the conduction bands with 
respect to the valence bands, to reproduce the correct band gap. The dielectric tensor in this 
approximation has been found to be within a few percent of the experimental data for more 
than a dozen semiconductors and insulators But, as emphasized by Dal Corso, Baroni, 
and Resta, there was no immediate justification within DFT of this successful procedure. 

In recent work by Godby, Sham and Schliiter [^], it was shown that the exchange- 
correlation potential near a semiconductor interface acquires a slow linear spatial variation, 
related to the discontinuity. A, that occurs on addition of an electron to the bulk semiconduc- 
tor. This variation corresponds to an ultra-non-local "vertex correction" in the Kohn-Sham 
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formulation of the dielectric response of the semiconductor, and depends crucially on the 
electron density at the interface [0. A connection between the scissor-operator correction 
and this linear behaviour of the exchange-correlation potential was also pointed out. 

In this paper we reexamine the response of a periodic solid to an electric field within 
DFT, and observe that the original proof of Hohenberg and Kohn does not apply to this 
case. We then prove the equivalent of the Hohenberg and Kohn theorem in perturbation 
theory, and show that knowledge of both the change in periodic density and the change in 
polarization are needed to recover the change of periodic potential and change of electric field 
that induced them. The exchange-correlation energy thus depends on the polarization, and 
this dependence, not taken into account previously, will generate an exchange-correlation 
electric field, identified with the slow linear variation of the exchange-correlation potential 
observed by Godby, Sham and Schliiter 0. We then explore a simple model that exhibits the 
crucial features of our theory, and explain the relationship to the scissor-correction approach 
to the dielectric tensor. 

In the DFT approach to the many-body problem, it is shown that the knowledge of the 
density n{r) of the ground-state of a system with Hamiltonian = T + Ve^e + (sum of 
the kinetic T, electron-electron interaction K-e, and one-body local potential v operators), 
uniquely determines the local potential v(y) of this Hamiltonian, up to a constant. 

Unfortunately, the straightforward application of an homogeneous electric field (linear 
potential), to a system with a periodic potential does not allow for a ground-state solution 
|[TT| : a translation against the direction of the field by a whole number of lattice constants 
would always lower the electronic energy. The impossibility of a ground-state in the presence 
of a finite electric field renders invalid the original proof of Density-Functional Theory 
for this case. 

The use of perturbation theory allows us to bypass this problem. We choose to work in 
the long wave method (an infinitesimal sinusoidal perturbation whose wavevector tends 
to zero). Equivalent results may be obtained using other techniques, mentioned in Ref. []T2| . 



Atomic units, for which the electronic charge is -1, are used. The change of total potential 



corresponding to an infinitesimal electric field 5£ is (written in one dimension for brevity) 

Sveir) = lim 5^.^^^^ = lim— (e'«^ - e-'^H. (1) 

We also allow for changes of potential that are periodic in space, with the same periodicity 
as the unperturbed system: bvcif) = Sv{G).e^'^^' with 6v{G) = {6v{—G))* where G is a 
non-zero vector of the reciprocal lattice. These changes of potential are obtained, in the 
long-wave method, from the Fourier components of the potential: 

Svcir) = \im{5v{G + q)e'^^+'^> + 5v{G-q)e'^^-''^'} (2) 

such that ^ = Sv{G + q) = Sv{G-q) . 

In response to these perturbations, at finite q, the system will develop changes in density 
described similarly by 6n{G±q). The long-wave part of this change in density, for g — * 0, 
will be III 



Snp{r) = lirn^P . gsin(gr), (3) 

where SV is the change of polarization for g = 0. 

The elaboration of a density functional theory for these perturbations must answer the 
following question: what quantities do we need in order to determine SS and the set of 6v{G) 
uniquely? 

Since we must stay within perturbation theory for treating electric fields, we now provide 
a perturbative analog of the first Hohenberg-Kohn theorem : the knowledge of the change 
in density everywhere uniquely determines the change in potential. This theorem could 
be proved by taking the infinitesimal limit of finite differences of the first Hohenberg-Kohn 
theorem [|1| , but this approach could not be followed for electric fields, since only infinitesimal 
electric fields are allowed. The following demonstration stays strictly within perturbation 
theory. 



We consider the Hylleraas minimum principle |T3|: a trial change in wavefunction 



gives an upper bound on the second-order change in energy 



5''E < {5(t>t\H - E\S(f>t) + {{S(l>t\Sv\(l>) + (c.c.)). (4) 

This principle is valid under the constraint + {(f)\S(f)t) — 0. The minimum is reached 

only for the 50 that is the response of the quantum-mechanical system to the change of 
potential 6v. When this change of potential is a 1-body local operator, Eq.(4) becomes 

6^E < {5(t)t\H - E\5(t)t) + j 6v{r)6nt{r)dr, (5) 

where the change in density Snt{r) is easily derived from the knowledge of the unperturbed 
wavef unction and the trial change in wavef unction. Now consider two changes in potentials 
5vi{r) and 5v2{r) such that 5vi{r) ^ 5v2{r) + constant. The Hylleraas minimum principle 
applied to the perturbation 6vi (r) gives 

{6(f)i\H - E\5(j)i) + j 5vi(r)5ni(r)dr < {5(t)2\H - E|(502) + j 6vi{r)6n2{r)dr, (6) 

while for the perturbation Sv2{r), a similar inequality, where 1 and 2 are interchanged, is 
obtained. Summing these two inequalities leads to 

< j {6vi{r) — 6v2{r)){6n2{r) — Sni{r))dr. (7) 

Setting Sni{r) = (5^2 (r) would lead to a contradiction, showing that two different changes 
in potential must induce two different changes in density. Thus, the knowledge of dn{r) 
everywhere uniquely defines the Sv{r) that induced it. 

The same line of argument can be used in the case of perturbations of periodic systems 
with finite wavevector q, as previously defined. All quantities have to be normalized to 
the unit cell volume. This normalization, and a Fourier transform, applied to the term 
/ 6v{r)6n{r)dr in Eq.(5), changes it into Q.ceu Y.ci^'v* {G+q)6n{G+q) + 6v*{G-q)6n{G-q)]. 
The limit g^O is now taken, for two different perturbations described by 5vi{G)} and 
{5£2,5v2{G)}. The G = term is isolated, and the long-wave values from Eq.(l) and (3) 
are used, such that : 

< ^ { {5£, - 5£2){SV2 - 5Vi) + E i^vUG) - 5v*2{G)){Sn2{G) - 5m{G)) }. (8) 

^ G^O 



If we now suppose 5Vi = SV2 and Sni{G) = Sn2{G), the expected contradiction is obtained. 
From this result we conclude that the change in potential and electric field can be deduced 
from the knowledge of the change in density and polarization that were induced by them. 
Note that the knowledge of the change of polarization is crucial, since it is the quantity 
conjugate to the change of electric field in Eq. (8) : if 6Vi were allowed to be different from 
SV2, Eq-(8) could be satisfied for some 6Si 7^ 682- The dependence on polarization is a 
remnant of the ultra-non-local dependence on the long-wave part of the change in density, 
Eq.(3), that vanishes for q = 0. 

Eq. (8) also shows that if no change in electric field is allowed, the knowledge of the 
polarization is not needed: a density-polarization functional theory is required for the specific 
case of the response to an electric field fl^. 

In the Kohn-Sham construction 0, an exchange-correlation energy functional of the 
density is introduced. Following the Kohn-Sham argument, now applied to the response to 
an electric field, this functional will depend on both the periodic density and the polariza- 
tion: E^d'i^G'i'P]- This functional is defined only for values of V in an infinitesimally small 
neighborhood of the zero-field polarization ]TB| . 

We show now that this polarization dependence of the exchange-correlation energy leads 
to the existence of a new "exchange-correlation electric field" . 

First, we mention that the treatment of the electric field as a perturbation makes use of 
the following link between the operator r and the derivative with respect to the wavevec- 
tor §: 

d 

S£ . Pcr\uj^k) = S£ . Pc'i'-Q^\uj,k) (9) 

where Pc is the projector on the conduction bands, and Uj^k is the periodic part of the Bloch 
function of band j with wavevector k. 



Also, the recent theory of polarization proposed by King-Smith and Vanderbilt [|1^] leads 
to the expression for the polarization [|T3| 
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fel^l^'t)^, (10) 

3 

where the sum is over the valence bands. 

When wavefunctions are varied in order to find the Kohn-Sham self-consistent equations, 
the exchange-correlation energy changes, owing to the implicit dependence of the density on 
wavefunctions, but also owing to the implicit dependence of the polarization on wavefunc- 
tions (see Eq.lO). The term Pc^Vxc\'Uj.k) that appears in the perturbative equivalent of the 
Kohn-Sham equation (the Sternheimer equation [^]) is replaced by 

d 

Pc5Vxc\Uj,k) + SS^c-Pci-Q^lUj^k) (11) 

where SS,, = ^SV + Eg ^Sug. 

The second term of Eq.(ll), missing in previous theories of the dielectric constant P|, 
is to be interpreted as an additional electric field (compare with Eq.(9)), created by the 
polarization dependence of the exchange-correlation energy. It vanishes in zero applied 
electric field. 

We will now analyze the model one-dimensional semiconductor used by Godby and Sham 
in Ref. [0, and show that it exhibits the main features of our theory: 1) the knowledge of 



the density change alone is not sufficient to deduce the periodic potential change and electric 
field change ; 2) an exchange-correlation electric field appears. 

In this model, the periodic solid plus electric field is treated by the long-wave method 
in a supercell consisting of N basic unit cells of length a. The external potential plus 
the Hartree potential Vext{x) + V}i{x) is taken to be the sum of two different contributions: 
Vo cos(^) +A sin(^). The first term has the periodicity of one unit cell, while the second, 
with the periodicity of the supercell, corresponds to a slowly varying potential of amplitude 
A. For an infinitely long supercell, the second term mimics the action of an electric field. 

The self-energy operator is taken to be the non-local potential: S(a;, x',a;) = 
_ x'\) where f{x) = — -Fo[l — cos{2nx /a)] is a negative function with the 
periodicity of one unit cell and g{y) is a normalized gaussian of width w = 2 a.u. 



We keep the same set of parameters as in Ref. |T0[- First, the many-body problem is 
solved by direct diagonalization of the equations containing the self-energy operator, using 
a plane- wave basis set and a Brillouin zone sampling at the F point only. From this result, 
we construct an exact density functional theory by determining the local potential Ves,i{x) 
which, when filled with non-interacting electrons (no self-energy operator), reproduces the 
same electron density as in the many-body case (see Fig. 4 of Ref. Il^l)- Standard iterative 



non-linear optimization techniques are used for that purpose. Independently, we also extract 
the (polarized) density of the single unit cell ±|a where the slowly varying potential is most 
linear and, again using optimization techniques, we reproduce this density with another 
potential, Ves,2{x), having period a. In this case, to retain the sampling of the Brillouin 
zone at the F point, we repeat the density periodically in the supercell. The only problem in 
reproducing the "target" density comes from the slight non-linearity of the slowly varying 
potential, and disappears progressively for larger supercells. Fig. 1 presents the results 
obtained with a 80a supercell and A= 0.95 eV, for which the r.m.s. difference between the 
target and computed density is 4 parts per thousand. 

This example illustrates that the action of a linear potential A(\4xt + Vr), superimposed 
on a periodic potential, in the case of an interacting-electron model system can be mimicked 
by two different potentials AV^g^i and AV^ff^2, acting upon the corresponding non- interacting 
electron system, in contrast to a naive application of DFT theorems. However, AVcs,2 is 
not able to reproduce the polarization (0.037 electrons), associated with the long- wavelength 
charge density. Moreover, A(V"ext + ^H) and the linear component AV^^l^^ of the Kohn-Sham 
potential AVefr.i differ by 15%, owing to the existence of the exchange-correlation electric 
field. 

We now discuss briefly the implications of our theory for practical calculations. LDA 
or GGA techniques can be used to approximate the periodic-density dependence of E^c- 
Unfortunately, since the homogeneous electron gas is metallic, we do not see any easy way 
to modify LDA or GGA in order to obtain a polarization-dependent exchange-correlation 
energy, that would generate the second term in Eq.(ll). Nevertheless, our theory allows 
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the scissor-operator correction to be understood as an approximate way of bypassing this 
dependence. 



Indeed, in Ref. ||I0|, it was shown that the effect of a slowly varying potential AV that 
induces transitions between valence and conduction bands with a gap Eg (see Eq.(lO) of 
Ref. [|I^) is approximately equivalent to the effect of a slowly varying potential AV + AVxc 
that induces transitions between valence and conduction bands with a gap E^^'^ (see Eq.(ll) 
of Ref. ||lO|). Hence, the exact density-polarization functional theory calculation in a field 
6S + 6£xc, with a gap Ej^^^, can be approximately replaced by a LDA calculation in a field 
6S, with a modified gap Eg-, this is the "LDA + scissor correction", which gives a dielectric 
constant for Si of 11.2 (Experiment 11.4-11.7, LDA 12.9-13.5, GGA 12.4-12.6) ||-|]. 

In conclusion, because the original Hohenberg-Kohn theorem is not valid for the case 
of the response to an homogeneous electric field, we have set up a density-polarization 
functional theory, of which an approximate form is the "LDA + scissor correction" . 
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FIGURE CAPTIONS 

Fig.l The unit cell of our model one-dimensional semiconductor, where the slowly varying 
applied potential A(\4xt + ^h), that changes the interacting electron density by An, 
is the most linear. Both AV^g i and AV"cfr.2, used in the non-interacting Kohn-Sham 
equations, yield the same An. AV^fj 2 is a periodic potential with no linear slope, while 
AKfr,i, whose linear part is AV^^f^, reproduces not only An but also the change of 
polarization due to A(V^xt + Vh). This illustrates the need for polarization-dependence 
in E;j.c- For clarity, the potential curves have been aligned so that they all start from 
zero. 
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